function [Hw, w, type, tao] = amplres(h)
    N=length(h);
    tao=(N-1)/2;
    L=floor((N-1)/2);
    n=1:L+1;
    w=[0:1000]*2*pi/1000;
    if all(abs(h(n)-h(N-n+1))<1e-10)
        Hw=2*h(n)*cos(((N+1)/2-n)'*w)-mod(N,2)*h(L+1);
        type=2-mod(N,2);
    elseif all(abs(h(n)+h(N-n+1))<1e-10)&&(h(L+1)*mod(N,2)==0)
        Hw=2*h(n)*sin(((N+1)/2-n)'*w);
        type=4-mod(N,2);
    else 
        error('错误，这不是线性相位FIR滤波器')
    end
    
end